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The star formation history of the LSB Galaxy UGC 5889 
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Abstract. We present HST photometry of the LSB galaxy UGC 5889 and derive its recent star formation history. 
In the last 200 Myr the star formation proceeded in modest bursts at a rate of the order of 10 _2 -10 -3 MQ/yr, 
with periods of extremely low SFR or even quiescence. The rate derived from the present study for the last 20 
Myr is in agreement with the H a emission from the galaxy. The presence of a consistent population older than 
200 Myr is suggested by the data. However, observational errors and completeness correction prevent any firm 
conclusion on the oldest age. The total mass of stars is of the order of 5.5 xlO 7 Mq. Even if the recent episodes 
of star formation have heated the gas and carved a hole in the disk, blow-away of the gas is unlikely to occur. 

Key words. Techniques: photometric; Stars: formation Galaxies: formation; Galaxies: stellar content 



1. Introduction 

It is well known that low surface brightness (LSB) galaxies 
contribute significantly to the total galaxy counts (see e.g. 
Marzke & da Costa, I1997f) and might even constitute the 
majority of galaxies (O'Neil & Bothun 2000). LSB galaxies 
show a large variety of properties, ranging from the blue 
to the red end of galaxy colours (O'Neil et al. 1997, 2000) 
and consisting of dwarfs as well as of massive systems 
(Sprayberry et al. 1995). Concerning their metal content, 
the vast majority of LSB galaxies show low metal con- 
tent (de Block &van der Hulst 1998), although some red 
LSB galaxies having solar metallicity have been recently 
detected (Bergmann et al. 2003). Among these diverse sys- 
tems, we focus on one of the most common types, the blue 
LSB (BLSB) galaxies (see the review by Impey & Bothun 
1997 and references therein), late- type disk galaxies hav- 
ing B surface brightness > 23 mag arcsec -2 where colours 
of (B-V)< 0.6 arc found (O'Neil et al. lTW7|) . These blue 
colours are puzzling primarily because of the restrictions 
they put on the star formation history of LSB galaxies. 
Even if they present low surface gas density, BLSB are 
among the most gas-rich galaxies at a given total lumi- 
nosity. The high gas fraction of BLSB galaxies, together 
with their low metal content, suggest that they are rela- 
tively unevolved objects. This can be interpreted in vari- 
ous ways. For instance BLSB galaxies can have very low 
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present and past star formation rates (SFR): in these sys- 
tems the main phase of star formation is still to occur. 
Alternatively, these galaxies can have a stellar popula- 
tion similar to those of high surface brightness spirals, but 
their luminosity can be dominated by a young population. 
Until now, the investigation of the stellar content of LSB 
galaxies has mainly been based on the discussion of the 
integrated colours. Many studies are consistent with a sce- 
nario where BLSB galaxies have on average a constant or 
increasing SFR with a ratio of young to old stars larger 
than found in high brightness spirals (van den Hoek et 
al. 2000, Knezek 1993, Bell & de Jong 1999). However, 
different opinions can be found in the literature, either 
suggesting that the SFR can be proceeding in frequent, 
small amplitude bursts (van den Hoek et al. 2000) or that 
the bulk of the population is older than 5-7 Gyr (Padoan 
et al. 1997). Finally, it is far from being understood how 
stars form in a low density, low metal content environment 
and what is the relation between SFR and other properties 
of LSB galaxies (i.e. the angular momentum) (see among 
others O'Neil et al. 2000, Boissier et al. 2003, Dalcanton, 
Spergel & Summers 1997). 

In this paper we focus on the study of a resolved LSB 
galaxy: UGC 5889 (NGC 3377A). This dwarf galaxy is a 
member of the Leo I group. Due to its similar redshift and 
resulting proximity - only 20 Kpc in projected distance 
(Sandage & Hoffman. 11991(1 . it is thought to be a compan- 
ion of the large elliptical NGC 3377. With M Hl /L B ~ 0.3 
Mq/ Lq, it is gas rich (Knezek et al. H999|l . It is classified 
as a low surface brightness dwarf having an average surface 
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brightness inside the Holmberg radius of B=23.97 mag 
arcsec ~ 2 . Moderately blue colours (U-B)= 0.06 ± 0.05, 
and (B-V) = 0.64±0.06 are given by Prugniel & Heraudeau 
<P§5|> . 

Base on its optical morphology, Knezek et al. (1999) 
suggested the presence of a spiral structure classifying 
the galaxy as SABm. This has recently been supported 
by high resolution spectroscopic H I-measurements of 
Huchtmeier et al. (2003 j wn0 nn d a double peaked pro- 
file for this galaxy which is typical for spirals. The HI 
intensity maps of Simpson & Gottesman (2000) show a 
ring-shaped distribution of the gas which is surrounding 
the bulk of the optical body and is correlated with the 
young star groups discovered by Vicari et al.(2002). 

Knezek et al. (|1999fl also measured the H a luminos- 
ity of UGC 5889 and found a relatively low value of 
Ljj a ~ 2.6 x 10 38 ergs s _1 , which is about ten times lower 
than the average H a luminosity in the sample of quies- 
cent dwarf spiral galaxies studied by van Zee et al. J1997). 
Thus, despite the dense environment where it is located, 
UGC 5889 appears to be a very quiet dwarf galaxy. Knezek 
et al. ((TTJ 99) discuss the possibility that it might represent 
a true transition object between late type dwarf galaxies 
and low SFR dwarf ellipticals as suggested by Sandage & 
Hoffman (1991). However they reject the hypothesis, since 
with its large gas reservoir UGC 5889 can continue to form 
stars at the present SFR for much longer than the Hubble 
time. 

The star formation history of UGC 5889 is unknown. 
It is still not clear whether the current SFR represents the 
mean level over the galaxy lifetime or, whether the star 
formation of this galaxy has occurred at a different rate 
in the past. This paper deals with the determination of 
the SFR and the SFR history from the analysis of colour- 
magnitude diagrams (CMDs). The plan of this paper is as 
follows: in Section|21 we discuss the observations, the data 
reduction and we present the observed CMD. In Section 
13 we derive the distance modulus, in section 0] the basic 
CMD features. In Section [S] we determine the star for- 
mation history of the galaxy and we discuss the feedback 
from star formation. In Section El we compare the SFR 
of UGC 5889 with other LSB galaxies. Conclusions are 
drawn in Section [7| 

2. Observation and reduction 

We use HST archive data of UGC 5889 taken in 1995 with 
WFPC2 at two different position angles and locations. In 
total, 21 partially overlapping F555W and four F814W 
frames are used including long and short exposed frames 
and yielding a total exposure time of 6000 s in F814W and 
25300 s in F555W. 

A composite image of the galaxy in the F555W filter 
is presented in Fig^ It shows the high spatial resolution, 
which allowed the photometric measurement of individual 
stars. 

The data have been reduced using the stand-alone ver- 
sion of DAOPHOT including ALLFRAME and following 



Table 1. The photometric errors as a function of the 
magnitude. 



F555W [mag] 


< 24 


25 


26 27 27.5 


o-phot [mag] 


< 0.1 


0.1 


0.2 0.3 0.45 


F814W [mag] 


< 24 


25 


26 


o-phot [mag] 


< 0.1 


0.2 


0.3 



the ALLFRAME cookbook (Turner EJ96) . To de-select 
suspected clusters and foreground galaxies we excluded 
the objects with FWHM larger than 2 pixels using the 
ALLFRAME parameter 'sharp'. In total, we derive pho- 
tometry for 4843 stars in the galaxy spanning the F814W 
magnitude range from 20 m to fainter then 27 m . 

The photometric calibration has been done following 
the prescriptions by Dolphin (|2000|l . We estimate the er- 
ror on the zero point as 0.1 mag. The completeness factors 
A are calculated as usual by means of artificial star exper- 
iments and are plotted in Fig. [21 The data are complete 
at 40% level for magnitudes brighter than F555W <27.5 
or F814W <26.5. Artificial star experiments are used as 
well to derive the errors on the magnitudes due both to 
the Poisson noise and crowding. These errors are listed 
in Table ^ Photometric errors due only to the Poisson 
statistics on the noise are estimated by ALLFRAME. As 
expected, they are smaller, being <0.1, 0.14, 0.22, 0.30 
for F555W magnitudes <25,26, 27,27.5, respectively. The 
analogous photometric errors on F814W magnitude are 
<0.1, 0.12, 0.24 for F814W = 24,25,26 respectively. 

The CMD of the whole galaxy is plotted in FigEI The 
blue plume is located at (F555W-F814W)< 0.5 and ex- 
tends from the limiting magnitude at F814W ~ 27 up to 
22. The red evolved stars extend from F814W ~ 27 to 
- 23 and have (F555W-F814W)>0.5. 

3. Distance modulus determination 

The I magnitude of the tip of the red giant branch (TRGB) 
has been shown to be an excellent distance indicator for 
nearby resolved galaxies (Lee et al.|1993). The position of 
the TRGB in UGC 5889 has been obtained dividing the 
observational F814W-(F555W-F814W) CMD into small 
intervals of width (0.15, 0.15) and then convolving the 
bi-dimcnsional histogram with an edge detector function, 
namely a gradient filter with an angle of 90 degrees, in- 
creasing the contrast at increasing magnitudes. The sharp 
cut-off at the end of the RGB is clearly detected by the 
gradient filter (see Fi^J) . We take the position of the main 
peak as our estimate of the tip magnitude and the half 
width at half maximum of the peak as the uncertainty 
on the determination. We find F814W=25.65±0.2 as the 
TRGB magnitude. Since all the existing calibrations of 
the TRGB magnitude are expressed in the Johnson pass- 
bands, we first convert the de-reddened HST F555, F814 
magnitudes to the V, I magnitudes. 
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Fig. 1. Composite image of the central area of the galaxy UGC 5889. 



We derive the foreground reddening from the maps 
of Schlegel et al. l(T55g|l . £? B -v = 0.034, resulting in 
^F555W-F8i4W — 0.06. This value is slightly lower than 
the -Bb-v = 0.06 estimated by Knezek et al. I|1999fl using 
the standard dust-to-gas ratio of Diplas & Savage l|1994[) 
and the H I measurements of Hartmann & Burton l|1997|l . 
The resulting difference in magnitude is within the pho- 
tometric errors and hence does not affect our further con- 
clusions. 

We therefore use £yb-v) = 0.034 corresponding to 
^(F555-F8i4) = 0.06 and A F 8i4 = 0.09. With these val- 
ues, the magnitude and colour at the TRGB transform 
into the Johnson I =25.47, (V-I)q=1.48. 

The distance modulus is derived from the Itrgb mag- 
nitude at the TRGB using the definition: 

(m — M)o = Itrgb — Mj^trgb 

The dependence of Mj RGB on the metal 
content [M/H] is given by the recent calibrations of 
Bellazzini et al.(2001, 2004): 

m trgb = .258[M/iJ] 2 + 0.676[M/H] - 3.626 

This holds true for populations in the age range from 
2 to 12 Gyr and -1.5 < [M/H] < -0.6. This method 



to derive the distance has been proved to be indepen- 
dent of the details of the star formation history (Barker 
et al. 2004). In the above metallicity range, 

M TRGB = 

-4.05 when [M/H] = -1.5 and Mj RGB = -3.91 when 
[M/H] = -0.6. Since the metal content of UGC 5889 
is likely inside the above range (see following Section), 
we are probably not seriously in error if we assume 
Mi : trgb = —3.98 ±0.07. This results in a distance mod- 
ulus of (m — M)o = 29.45 ± 0.2, where the quoted errors 
take into account both the uncertainties on the absolute 
magnitude and on the observational magnitude of the tip. 
This value is in reasonable agreement with previous deter- 
minations. For the companion NGC 3377, Graham et al. 
(1997) derive a distance modulus 30.01 ±0.15. A direct de- 
termination for UGC 5889 has been made by Makarova & 
Karachentsev II1998|I using the magnitude of the brightest 
blue stars, and they find (m-M) =29.8. 

4. CMD discussion 

To derive information about the SFR, we compare the 
CMD with isochrones of different metallicity. Fig|S] shows 
the CMD with the isochrones of Z = 0.008 which repro- 
duce the main features of the CMD such as the slope of 
the giant branch and the location of those stars brighter 
than F814W- 26 and bluer than (F555W-F814W)~ 0.5 
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Fig. 2. Completeness factors A for F555W and F814W 
magnitudes 
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Fig. 3. CMD of UGC 5889 (whole field) 



Fig. 4. RGB tip detection: the gradient filter response 
to the luminosity function is plotted as a function of the 
magnitude(see text for more detail). The position of the 
main peak is taken as indication of the RGB tip (solid 
line). The dotted lines show the half width at half max- 
imum of the peak taken as estimate of the uncertainties 
on the tip detection 



about 200 Myr. At this age, main sequence stars and/or 
He-burning loop stars are brighter than the completeness 
limit. 

Concerning the presence of an older population, hints 
can be derived from the CMD. The large clump redder 
than (F555-F814)~ 1.2 extending down to the photometry 
limit of F814W ~ 27 is consistent with the presence of 
AGB and He-burning stars of Z = 0.008 in the age range 
~ 600-1000 Myr. Comparing the observational CMD with 
isochrones as metal-poor as Z = 0.0001, it is evident that 
this metallicity cannot reproduce the AGB/RGB colour 
of stars younger than 1 Gyr (see Fig.EJ). However, the 
existence of a metal poor population ( Z < 0.0001) in 
the age range 1-12 Gyr cannot be ruled out. Stars of this 
age in the RGB/ AGB phase are expected to appear at 
magnitudes fainter than F814W ~ 23.5 and at a colour 
(F555-F814)> 1.2. 



(blue plume). This blue plume is formed by a mixture of 
main sequence and blue loop stars which crowd the same 
CMD region. 

A solar metallicity is only marginally consistent, pre- 
dicting a colour of the core He-burning stars redder than 
observed (see Fig|HJ|. Core He-burning stars as old as 200 
Myr are crowding the clump of the red stars brighter 
than F814W - 27 and having (F555-F814)< 1.2. We ex- 
pect to have a good time resolution for ages younger than 



5. Star formation rate 

In order to infer the SFR of this galaxy, theoretical CMDs 
in different age ranges are simulated. These include the 
spread due to observational photometric errors. For each 
age bin, 1000 stars were generated. The generation of the 
synthetic populations makes use of the set of stellar tracks 
by Girardi et al. iTH^ l27TO]l . The initial mass function 
(IMF) of Kroupa (2001, 2002) is assumed. This IMF is a 
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Fig. 5. CMD of UGC 5889 with isochrones of metallicity 
Z = 0.008. For clarity, not all the stars in the field of UGC 
5889 are plotted. 



Fig. 7. CMD of UGC 5889 with isochrones of metallicity 
Z = 0.0001. 
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Fig. 6. CMD of UGC 5889 with isochrones of metallicity 
Z = 0.02. 



power-law function with a slope a — 2.3 for stellar masses 
m > 0.5 Mq, while a — 1.3 in the mass range 0.08-0.5 
Mq (when the standard Salpeter value is 2.35). The con- 
stancy of the IMF slope in different environments is still 
a matter of discussion, although a number of recent pa- 
pers proposed the idea of a universal IMF (see e.g. Kroupa 
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Fig. 8. The subdivision of the CMD into bins. 



2002, Wyse et al. 2002, Chabrier 2003, Weidner & Kroupa 
2004). Recently, it has been suggested that a might de- 
pend on the galaxy type. For the IMF in LSB galaxies, it 
has been proposed that the effect of the low metal content 
can result in an inefficient SF taking place mainly outside 
the giant molecular clouds, thus yielding an IMF deficient 
in massive stars. As a consequence the IMF slope might 
be as steep as a ~ 3.8 (Schombert et al. 1990). However, 
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in this paper a more conservative approach is used, since 
a variation of the IMF slope in LSB galaxies is not at 
present completely justifiable (Lee et al. 2004). The com- 
pleteness of the data is taken into account by dividing the 
simulated CMD in magnitude-colour bins and then sub- 
tracting from each bin having stars, (l-A)N t ^, where 
A is the smallest of the F555W and F814W completeness 
factors as given in FigO 

Finally, the SFR is derived by means of a downhill sim- 
plex method (Harris & Zaritsky 2001, Rizzi et al. 2002), 
minimising the \ 2 function in a parameter-space having N 
dimensions. At each step the local \ 2 gradient is derived 
and a step in the direction of the gradient is taken, until 
a minimum is found. In the following, the observational 
CMD is divided into bins. Recent work concerning the de- 
termination of the SFR from the CMDs has pointed out 
the importance of using binning that takes into account 
the various stellar evolutionary phases, as well as the un- 
certainties on the stellar models (Rizzi et al. 2001). For 
this reason, while a coarser magnitude bin distribution is 
used, only two bins in colour are considered, namely from 
-0.5 to 0.5, corresponding to blue plume stars, and from 
0.5 to 3, where the evolved stars are located. This avoids 
that the uncertainties on both the observational errors 
and the theoretical models (i.e. on bolometric corrections, 
RGB and AGB location, extension of the core He-burning 
loop) result in spurious solutions. Additionally the mag- 
nitude bins are taken in the blue part of the CMD at 
constant F814W magnitude, while in the red part they 
are at constant F555W magnitude. The reason for this is 
that the completeness limit in the red part of the CMD 
is mainly dependent on the completeness of the F555W 
magnitude. This is illustrated in Fig. [8] To prevent set- 
tling on local rather than global minima, the simplex is 
first started from a random position, then when a possi- 
ble solution is obtained, it is re-started from a position 
very close to it. Finally, when a minimum is found, 30000 
random directions are searched for a new minimum. 

The first guess solution is obtained by comparing the 
observational CMD with isochrones of different ages and 
metallicities (see previous section). We use N=14 stellar 
populations, whose ages are listed in Tabled The chemical 
enrichment law is assumed as a pre-defined input parame- 
ter, derived comparing data and isochrones. Following the 
discussion of the previous section for ages younger than 20 
Myr, we assume that the stars are generated with a metal 
content randomly chosen in the range Z=0. 004-0. 01. SFR 
episodes in the age range 20 Myr-200 Myr are assumed to 
be produced by stars having Z=0. 001-0. 01. An older pop- 
ulation in the age range 0.2-1 Gyr with the same metal 
content is included. A population of metal-poor stars with 
Z=0. 0003-0. 004 and as old as 1-12 Gyr is also taken into 
account, since due to the low metal content, their magni- 
tudes will be brighter than the photometry limit of F555W 
~ 27.5. Metal-rich stars in the same age range are not 
considered since they are expected to be too faint. As dis- 
cussed in Section 0] we are aware that we have poor age 
and SF resolution for populations older than 200 Myr. 
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Fig. 9. SFR of UGC 5889 as a function of look-back time. 
Note that due to the low age resolution beyond lGyr, we 
only give the average SFR in the lower plot (see text and 
table EJ). 

Table 2. The solution parameter 





tin 


6 fin 


SFR 


SFRi ow 






Gyr 


Gyr 


lO- 3 M /y 


lO" 3 M /y 


1O- 3 M /j/ 


SFR1 


0.006 


0.007 


86.9 


85.3 


90.9 


SFR2 


0.010 


0.015 


1.6 


1.6 


2.1 


SFR3 


0.015 


0.020 


3.7 


3.6 


4.2 


SFR4 


0.020 


0.030 


0.0 


0.0 


0.6 


SFR5 


0.030 


0.040 


0.17 


0.17 


0.8 


SFR6 


0.040 


0.050 


2.0 


1.7 


2.8 


SFR7 


0.050 


0.060 


0.0 


0.0 


0.9 


SFR8 


0.060 


0.070 


0.0 


0.0 


1.0 


SFR9 


0.070 


0.080 


6.9 


6.7 


7.9 


SFR10 


0.080 


0.090 


0.0 


0.0 


l.l 


SFR11 


0.090 


0.100 


0.0 


0.0 


l.l 


SFR12 


0.100 


0.200 


7.0 


6.9 


7.3 


SFR13 


0.200 


1.000 


1.8 


1.7 


2.2 


SFR14 


1.000 


12.000 


7.5 


7.3 


7.6 



Therefore, the old age bin is included only to get an indi- 
cation of the presence of an old population. 

In Fig|5|the best solution is plotted. The correspond- 
ing best fit SFRs are given in Table [21 We derive a \ 2 
value of 1.7. In order to better understand the quality of 
the fit, we compare the values of \ 2 from two simulations 
with the same observational errors and input parameters, 
but different random number generators. The resulting 
X 2 for the second solution is 1.5 and supports the orig- 
inal finding. The range of SFR corresponding to 68% of 
the confidence level is given in Table This formal error 
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is obtained by exploring the parameter space by varying 
all the amplitudes simultaneously in random directions, 
stepping away from the minimum until the \ 2 value in- 
dicates that we have reached the 68% confidence level. 
The number of random directions to be checked has been 
tested on the simulations. About 50000 random directions 
are explored and are found to be sufficient to derive the 
confidence level. 

The recent star formation appears to have proceeded 
in modest episodes at a rate of the order of 10~ 2 MQ/yr, 
with periods of lower star formation or even quiescence. 
A very short burst took place 6 Myr ago and lasted about 
1 Myr. Any other solution involving different durations, in- 
tensities, or age scales for star formation episodes younger 
than 10 Myr yields a worse fit. 

To assess the reliability of the solution at old ages, 
the downhill simplex method is applied to a different 
sets of initial ages, excluding in turn the oldest popu- 
lations, in the age range 0.2 - 1 Gyr and 1-12 Gyr, re- 
spectively. The simplex converges at high values of x 2 > 
4.6 and 2.9 respectively. To better verify the presence 
of an old population, we consider the CMD stars redder 
than (F555W-F814W)- 0.5 and fainter than F555W~ 26. 
This region of the CMD is populated by stars older than 
0.09 Gyr. We analyse the colour distribution in two mag- 
nitude bins, namely 26 < F555W< 27, 27 < F555W< 
27. 5, corresponding to the two faintest magnitude bins in 
the red part of the CMD in FigO We make use of the 
last 4 populations of Table [2] Figs EH and ^] present the 
best fit solution together with the data in two cases, a) in- 
cluding the population older than 1 Gyr, and b) without 
it. In the case b) it is not possible to fit all the magni- 
tude bins with the same combination of populations. For 
instance in Fig. ^] the SFR fitting the bin 27-27.5 can- 
not reproduce the data in the brighter magnitude bins. 
This is evident as well from Figll2l where the confidence 
interval for the SFRs of populations 12 and 13 are plotted 
in case b). No common solution is found when the pop- 
ulation older than 1 Gyr is not included. However, tak- 
ing into account the presence of observational errors and 
completeness correction, we prefer to conclude than the 
resolution of the data is sufficient to assess the presence of 
a consistent population older than 200 Myr, but finer age 
division cannot be reliably derived. Finally, the stability 
of the solution against different sub-divisions of the CMD 
in cells is tested: slightly larger uncertainties result when 
a constant bin of 0.5 mag is used. 

Fig El presents the simulated CMD for the best solu- 
tion. The simulated CMD is similar to the observational 
one. The total mass going into stars is about 5.5 xlO 7 
Mq. Since the total mass of the gas at present is ob- 
servationally estimated to be 1.6 xlO 7 Mq (Simpson & 
Gottesman |2 Q0 _ Qf . the present day gas fraction comes out 
to be fg=M gas /(M gas + M s t a rs) = 0.2. An independent 
check of this value can be derived if i g is expressed in 
terms of observable quantities as: 

fg = (1 + M stars /L B x Lb/tiMhi)- 1 

Bell et al. (2003) derive for a Salpeter IMF: 



,= uu : 26< F555W < 37 

150 ' T I . 

100 r-*- - L % \ ' 1 

"^>--f-4-"°"T-T"r- it "-i~" ~ l t ^"h ^ 

300 - 27< F656W < 27.5 



200 - 



100 

0.5 ~ 1 1.5 2 2.5 

(F555W F814W) 



Fig. 10. The observational colour distribution (heavy 
solid line) compared with the simulations (thin solid line) 
in two magnitude bins when the population 13 older than 1 
Gyr is taken into account. Open squares indicate the sim- 
ulated colour distribution in the age range 0.09-0.1 Gyr, 
crosses are the analogous values in the age range 0.1-0.2 
Gyr. Filled squares show the analogous values in the age 
interval 0.2-1.0 Gyr. Open triangles represent the simu- 
lated distribution for ages older than 1 Gyr. 

log(M stars /L B ) = 1.737(5 - V) - 0.789 
Using (B-V)=0.64, M HI /L B = 0.3, and rj = 1.4 we 
find f g = 0.17 in agreement with the determination ob- 
tained from the SF. 

5.1. Star formation rate from H a luminosity 

The total luminosity ~LH a emitted in H Q can be used to 
estimate the recent SFR. Following Kennicutt et al. Q1994[l 
the SFR can be obtained from 

SFR = 1.26 x 10~ 41 x L Ha /l0 3S M^yr' 1 . 

Adopting Ljj a ~ 2.6 x 10 38 erg s" 1 as derived for 
UGC 5889 by°Knezek et al. we find SFR ~ 4.0 x 

10~ 3 M Q yr _1 . This SFR should be regarded as the av- 
erage value over the lifetime of the stars emitting in H a . 
Assuming that the majority of the flux in H Q is due to 
stars more massive than 10 Mq, this gives an estimate 
of the SFR over the lifetime of a 10 M Q star. This life- 
time can be derived from Girardi et al. (2000) tracks for 
Z = 0.008 as 20 Myr. This SFR can be directly compared 
with the rate derived from the star-counts. Using the coef- 
ficients of Tabled we derive a SFR of 5.6 x 10~ 3 Moyr" 1 
averaged over the past 20 Myr of the galaxy lifetime, in 
good agreement with the ~Ljj a determination. 
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Fig. 11. The observational colour distribution (heavy Fig. 12. The confidence interval for the SFR of popula- 



solid line) compared with the simulations (thin solid line) 
in two magnitude bins when the population 13 older than 
1 Gyr is not included. Open squares indicate the simu- 
lated colour distribution in the age range 0.09-0.1 Gyr, 
crosses are the analogous values in the age range 0.1-0.2 
Gyr. Filled squares show the analogous values in the age 
interval 0.2-1.0 Gyr. 



tion 12 and 13 derived from the red stars in the faintest 
magnitude bins, namely 26 < F555W< 27 (1), 27 < 
F555W< 27.5 (2), when no population older than 1 Gyr 
is included in the data (see text for details). 



20 



5.2. The feedback from star formation 

The HI distribution of UGC 5889 shows a concentration 
of the gas in a ring surrounding the bulk of the optical 
body (Simpson & Gottesman 2000). The velocity field re- 
veals a rotation of the entire ring-shaped distribution and 
no sign of expansion of the ring is found. This kind of dis- 
turbed interstellar medium is not unusual in dwarf galax- 
ies (Puche & Westphal 1993, Young & Lo 1996, Simpson 
& Gottesman 2000). The hypothesis naturally arises that 
this dramatic morphology might be the result of a vig- 
orous star formation episode in a galaxy with a shallow 
potential well. In fact the energy injected in the gas by 
SNs and stellar winds can trigger the formation of super- 
bubbles that can easily expand in a low density medium. 
The expansion sweeps up the gas, forms a cavity at the 
centre, and leaves behind a system where the majority of 
the gas is located in the outskirts. This gas can be lost by 
the galaxy or can cool and rain back on to the galaxy after 
a cooling time (McLowfe Ferrara 1999). In this section, we 
discuss whether the energy injected in the gas by the SN 
explosions we estimate on the basis of the derived SFR 
can be sufficient to explain the morphology of the HI in 
UGC 5889. Finally, we discuss what might be the fate of 
this gas, i.e. whether a blow-out or a blow-away can take 
place. Following McLow& Ferrara (1999) and Ferrara & 
Tolstoy (2000) we define a blow-out as a process where 



Simulated CMD 



22 



24 



26 




2 4 
(F555W FB14W) 



Fig. 13. The simulated CMD using the SFR of Table El 



the central SN explosions carve a hole in the gas disk, ac- 
celerating a fraction of it, while the remaining gas stays 
almost unperturbed. In a blow-away nearly all the gas is 
accelerated above the escape velocity and is lost from the 
galaxy. 
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To discuss the possibility of blow-out, let us assume 
that, perpendicular to the galactic plane, the unperturbed 
gas density distribution in the galaxy is exponential p(z) = 
Poexp(— z/H), where z is the vertical coordinate, H is the 
galactic exponential scale height. The velocity of the shock 
wave produced by the SN explosions can be expressed as 
v ~ (p/p)' 1 / 2 ), where p is the pressure. The velocity v 
of the shock produced in a stratified medium decreases 
down to a minimum occurring at z=3H before being re- 
accelerated in case of a blow-out. A blow-out takes place 
when the velocity exceeds the escape velocity of the galaxy 
at z=3H, where the shock wave is accelerated to infinity. 
LSB galaxies are believed to be dark matter dominated 
(McGaugh & de Block 1998). In order to calculate the 
escape velocity, the galaxy is modelled including the pres- 
ence of a dark matter halo as in Persic & Salucci (1996). 
The dark to visible mass ratio <p can be expressed as a 
function of the visible galaxy mass M g ^ in 10 7 Mq, as 

<j) ~ 34.7M~ 7 - 29 . 

Assuming as M g j = 19 (Simpson & Gottesman 2000), 
we find 4> ~ 14. This value is in reasonable agreement with 
the ratios in the range 10-30 derived for LSB galaxies by 
McGaugh & Blok (1998). The blow-out condition can be 
expressed as (McLow & Ferrara 1999): 

L 3S > 1.2 x 10- 2 M g>7 c 2 w h4>/(ct> + 1) 

where L38 is the mechanical luminosity of the SN ex- 
plosions i.e. the energy produced by SN explosions divided 
by the average lifetime of the less massive star m; ou) ex- 
ploding as a SN. If we assume mi ow ~ IOMq then its life- 
time is 20 Myr. L 3 s is expressed in 10 38 erg sec -1 when 
cio — 1 is the sound speed of the interstellar medium in 10 
km sec -1 including a turbulence contribution, and h is the 
scaled Hubble constant, assumed to be 0.65. Comparing 
the HI location with the optical image, it is clear that the 
youngest stars are located in the outskirts of the gas dis- 
tribution. Hence, the youngest star formation episode at 
5 Myr cannot have induced the SN explosions responsible 
for the central gas depletion. For this reason, we calculate 
L38 using the mean SFR derived from our model in the 
last 20 Myr. Substituting in the blow-out condition we 
obtain: 

L 38 ~ 6 >> 0.14 

This means that the blow-out condition holds true in- 
side the gas: a blow-out is then responsible for the peculiar 
morphology of the ISM in UGC 5889. By definition, the 
blow-out involves only a fraction of the mass of the gas, 
that which is inside the cavities created by the SN explo- 
sions. A more disruptive phenomenon is the blow-away in 
which the gas content of the parent galaxy is completely 
lost. The blow-away can take place when the momentum 
of the gas shell at the radius r s is larger than the momen- 
tum necessary to accelerate the gas outside r s at a velocity 
larger than the escape velocity. This condition can be ex- 
pressed following McLow & Ferrara (1999) as: 



L 38 > 8 x 10- 2 M^/a, o )V 10 ft 

where uj = 3 is a scale factor. Under reasonable as- 
sumptions concerning the sound speed c%o — 1.5 — 0.5, 
this condition is not verified even if <p is assumed to vary 
in the range 5-30. The ultimate fate of the gas seems to 
be that it will not leave its parent galaxy. This result is in 
full agreement with the low gas velocity derived observa- 
tionally (Simpson & Gottesman 2000). 

6. Comparison with other LSB galaxies 

Bursty star formation history in LSB galaxies has been 
suggested in the literature. We quote, from amongst oth- 
ers, Boissier et al. 1)2003(1 who find that it is necessary to 
include star formation events such as bursts in the mod- 
els of LSB galaxy evolution to account for the observed 
scatter in the Tully-Fisher and total mass of the gas-to- 
luminosity relation. A burst scenario is proposed by O'Neil 
et al. 1)1998)1 to explain both the blue colour as well as the 
disrupted morphology of most blue LSB galaxies. 

The dwarf irregular galaxy IC 1613 has similar colours 
(B-V = 0.71) and its SF has been studied in great detail. 
Cole et al. (1999) find a nearly constant SFR at a low 
level of 1.6 x 10" 3 M Q yr _1 kpc" 2 over the past 1-300 Myr. 
However, due to the large angular size of IC 1613, they 
only study the very central part of it, which might not be 
representative of the SF over the whole galaxy. Note for 
example that in UGC 5889, the recent SF is concentrated 
on the outer edge, while the centre is inhabited by the 
older population. Skillman et al. I|2003|l compared it with 
a field situated at the outer edge of IC 1613, but studied 
the intermediate-age and ancient SFR, starting at an age 
of 100 Myr. A comparison with the CMD of Cole et al. 
yielded no significant difference. Although their data are 
consistent with constant SFR, they find small bursts in the 
SFR and a strong enhancement of the SFR at intermediate 
ages. However, no evidence is found for a strong, bursty 
SFR in the recent history. In fact, the CMDs presented 
by Skillman et al. as well as the one by Cole et al. lack 
the bright blue stars that are present in UGC 5889 and 
can only be explained by recent SF. Therefore, although 
both galaxies are classified as blue LSB galaxies they differ 
in their SFR in the sense that IC 1613 is dominated by 
SF at intermediate age, while UGC 5889 is dominated by 
more recent SF. This behaviour is also reflected by the 
total mass-over-luminosity ratio, which is slightly lower 
for UGC 5889. 

7. Conclusions 

In this paper we have discussed the star formation his- 
tory of the LSB galaxy UGC 5889 on the basis of HST 
archive data. Our main conclusions can be summarised as 
the following: 

1) the RGB tip is detected at F814W=25.65±0.2 im- 
plying a distance modulus of (to — M) = 29.45 ± 0.2. 
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2) the recent star formation is proceeding in modest 
bursts at a rate of the order of 10~ 2 - 10~ 3 M /yr. This 
rate is lower than the 5-10 Mq/jt derived for high sur- 
face brightness disk galaxies (Kennicutt 1992), but only 
slightly larger or comparable to the 1CT 3 MQ/yr observed 
in dwarf irregulars (Hunter & Gallagher 1986). The SFR 
derived in UGC 5889 from H Q emission is in agreement 
with the results obtained from CMD analysis. 

3) The existence of a consistent population older than 
200 Myr and possibly older than 1 Gyr is suggested by our 
analysis. However, this result might depend on the com- 
pleteness correction and observational errors and needs to 
be verified with deeper photometry. 

4) The present-day total mass in stars is of the order 
of 8.3 xlO 7 Mq. Since the total mass of the gas has been 
estimated as 1.6 xlO 7 Mq, the gas fraction is f g =0.19. 

5) The morphology of the gas inside this galaxy re- 
veals a large cavity, including the main optical body of the 
galaxy. The gas is mainly located in the external regions. 
We have calculated the SFR feed-back on the interstellar 
gas. We find that this impressive morphology is consistent 
with the hypothesis that SN explosions have produced a 
blow-out of the gas from the internal regions towards the 
external ones. On the basis of our calculations, it seems 
however quite unlikely that the final fate of this gas is to 
leave the galaxy. 
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